MALACOLOGIA, 2013, 56(1-2): 231—252 


MOLECULAR PHYLOGENY AND BIOGEOGRAPHY OF A HIGH MOUNTAIN 
BIVALVE FAUNA: THE SPHAERIIDAE OF THE TIBETAN PLATEAU 


Catharina Clewing1*, Ulrich Bóssneck?, Parm V. von Oheimb1 & Christian Albrecht! 


ABSTRACT 


Sphaeriids represent a highly diverse family of small bivalves inhabiting freshwater habitats 
worldwide. They are recognized for their high dispersal capacity as well as for a high capa- 
bility to adapt to such extreme environments as high mountain regions over 4,000 m above 
sea level (a.s.l.). The Tibetan Plateau, located in central Asia, is such a high elevation area. 
Here, we present the first broad spatial sampling of sphaeriids from the Tibetan Plateau, 
which allowed us to describe the Recent fauna as well as their biogeographical background. 
Our study revealed six sphaeriid species inhabiting various habitats on the Tibetan Plateau: 
Pisidium nitidum, P. stewarti, P. subtruncatum, P. zugmayeri, Pisidium sp. | and Musculium 
kashmirensis. Four species are newly recorded for this region. Based on phylogenetic analy- 
ses using molecular information from three different genes (COI, 16S and 28S) we conclude 
that the extant diversity of Tibetan Plateau sphaeriids might have been caused by multiple 
colonization events. Furthermore, the zoogeographic affinities of Tibetan Plateau sphaeriids 
are mostly related to the Palearctic. Similar patterns have been observed in previous studies 
of other freshwater molluscan taxa in the area. 
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INTRODUCTION 


The family Sphaeriidae comprises small to 
minute freshwater bivalves with shell lengths 
ranging from 2.5 to 25 mm (Korniushin, 2006). 
They are the second most diverse family of 
freshwater bivalves after the Unionidae, with 
an estimated number of 227 species (Graf, 
2013) in five genera: Eupera Bourguignat, 
1854; Byssanodonta d'Orbigny, 1846; Spha- 
erium Scopoli, 1777; Musculium Link, 1807; 
and Pisidium C. Pfeiffer, 1821 (Dreher Mansur 
& Meier-Brook, 2000; Korniushin & Glaubrecht, 
2002). However, phylogenetic relationships 
among and within the sphaeriid taxa are still 
largely unresolved (Lee & O Foighil, 2003; 
Korniushin & Glaubrecht, 2006, and references 
therein; Schultheiss et al., 2008), which results 
in conflicting classifications especially on the 
generic level (Graf et al., 2013). 

ophaeriids are currently distributed worldwide 
(except for Antarctica; Bogan, 2008) and are 
known for their high capability to adapt to di- 
verse habitats and extreme environments, such 
as high mountain regions > 4,000 m above sea 


level (a.s.l.) (Preston, 1909; Prashad, 1933; 
Bóssneck, 2012). Furthermore, sphaeriid spe- 
cies may have a large distribution range (i.e., 
transcontinental distribution), but can also be 
restricted to small ranges (Korniushin, 2006). 

These characteristics as well as their high 
passive dispersal capacity (Kappes & Haase, 
2012) make the Sphaeriidae suitable model 
organisms for phylogeographic studies, par- 
ticularly in species-poor high-elevation regions 
such as the Tibetan Plateau. 

The Tibetan Plateau of central Asia is the 
largest and highest plateau on earth (mean el- 
evation of ~5,000 m a.s.l.; Royden et al., 2008) 
and surrounded by high mountain ranges: Qilian 
and Altun Mountains in the north, Himalaya and 
Karakorum Mountains in the south and west, 
and Longmen Mountains in the east (Chang et 
al., 1989). Approximately 1,600 freshwater and 
saline lakes greater than 1 km? are scattered 
throughout the plateau (Zheng, 1997). These 
lakes are typically oligotrophic (Xing et al., 
2009), and many of them are part of endorheic 
drainage basins (Revenga et al., 1998). Besides 
the lakes, numerous smaller lentic water bodies, 
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such as ponds and wetlands, are present (An | 


et al., 2009). The plateau with its large glaciers 
("third pole"; Xu et al., 2009) is of great impor- 
tance for regional freshwater systems, because 
several major rivers of South and East Asia 
(such as the Yangtze, Mekong, Brahmaputra, 
and Indus) have their headwaters in this area. 

The faunal components of the Tibetan Pla- 
teau are mostly thought to be of Palearctic 
and Oriental origin (Yang et al., 2009), and the 
freshwater mollusk fauna seems to be mainly 
characterized by Palearctic species (Oheimb et 
al., 2011, 2013). Thereby, the Himalayan moun- 
tain range forms the southern border to the 
Oriental region (Cox, 2001; Morrone, 2002) and 
is assumed to be a major dispersal barrier for 
freshwater taxa (Bánárescu, 1992; Abell et al., 
2008; Oheimb et al., 2013). In addition, extant 
patterns of diversity and distribution of species 
throughout the Tibetan Plateau are mainly in- 
fluenced by extreme geological changes during 
the plateau's uplift, as well as climate fluctua- 
tions especially during the Pleistocene (e.g., 
Last Glacial Maximum, ca. 20,000 years ago; 


Yang et al., 2009). Although the Tibetan Plateau 
was probably not completely covered by an 
ice sheet throughout the Pleistocene (Zheng & 
Rutter, 1998; Shi, 2002), the unstable climatic 
conditions resulted in massive environmental 
changes over time, such as lake-level fluctua- 
tions (Lehmkuhl & Haselein, 2000; Dietze et 
al., 2010). These changing environments might 
have triggered diversification processes in 
Tibetan Plateau freshwater organisms, though 
the regional freshwater fauna and particularly 
the mollusk fauna remains species-poor and 
scarcely studied. A total of three gastropod 
genera (Radix Montfort, 1810; Gyraulus 
Charpentier, 1837; Valvata O. F. Müller, 1773) 
and two genera of bivalves (Musculium, Fig. 
1; Pisidium, Figs. 2-6) were reported for the 
plateau (Prashad, 1933; Mischke et al., 2010; 
Bóssneck, 2012; Oheimb et al., 2011). 

Two recent studies focused on the biogeo- 
graphical history of pulmonate freshwater 
gastropods (the lymnaeid genus Radix and 
the planorbid genus Gyraulus) on the Tibetan 
Plateau (Oheimb et al., 2011, 2013). 


FIGS. 1-6. Sphaeriid species sampled across the Tibetan Plateau. FIG. 1: M. kashmirensis (I), location 
CN19; FIG. 2: P. stewarti (VI), CNO7; FIG. 3 Pisidium sp. | (II), CNO7; FIG. 4: P. nitidum (MI), CN21; 
FIG. 5: P. subtruncatum (IV), CN14; FIG. 6: P. zugmayeri (V), CN22. Roman numerals refer to the 
respective phylogenetic lineages, see Fig. 12. Scale bar = 0.5 mm. 
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Although sphaeriid clams often occur in 
sympatry with Radix spp. and Gyraulus spp., 
major biological and ecological differences 
exist, potentially affecting phylogeographical 
patterns. First, the clams are soft-substratum- 
dwelling filter-feeders with a potential to 
influence ecosystem processes in freshwater 
systems (Vaughn & Hakenkamp, 2001; Dil- 
lon, 2004). Secondly, sphaeriid species show 
the reproductive strategy of euviviparity, that 
is, eggs, embryos, and juveniles are bred in 
specialized structures of the parental gills 
(Korniushin & Glaubrecht, 2006), providing a 
comparatively stable environment for the off- 
spring. Thirdly, the (passive) dispersal capac- 
ily of sphaeriids is presumed to be very high 
(Mackie, 1979; Schultheiss et al., 2008), thus 
playing a major role in colonizing new habitats 
(Zelaya & Marinone, 2012). In fact, there are 
several mechanisms enabling these small 
bivalves to disperse between waterbodies in 
close proximity, as well as over long distances 
(Kappes & Haase, 2012). Vector-mediated 
passive dispersal seems to be of particular 
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importance. For example, clams were found be- 

ing attached to the limbs of freshwater insects 

and amphibians (aquatic vectors; Rees, 1965; 

Bilton, 2001) and to the feathers and legs of dif- 

ferent waterbirds, as well as to the extremities 

of flying insects (aerial vectors; Wesselingh et 
al., 1999; Figuerola & Green, 2002; Zelaya & 

Marinone, 2012). Furthermore, Mackie (1979) 

reported that juvenile sphaeriids are able to 

survive the gut passage in ducks. 

Considering these particular characteristics of 
Sphaeriidae (in contrast to other Tibetan Plateau 
freshwater mollusks) as well as the special en- 
vironmental conditions of the high-elevation Ti- 
betan Plateau, this study has two major aims: 

(1) Assessment of the Tibetan Plateau Sphaerii- 
dae fauna (including data about diversity and 
ecology) based on extensive sampling. 

(Il) Reconstruction of the biogeographical back- 
ground of the Sphaeriidae species inhabit- 
ing the plateau (e.g., number of colonization 
events, zoogeographic affinities) based on 
the phylogenetic information inferred from 
three different DNA markers. 


FIG. 7. Sampling sites of Sphaeriidae across the Tibetan Plateau and Nepal. Locality information is 


given in the Supplementary material (Appendix 1). 
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MATERIAL AND METHODS 
Sampling 


During field work conducted between 2009 
and 2012, specimens of the family Sphaeriidae 
were found at 27 locations throughout the Tibet- 
an Plateau (Fig. 7). Specimens were collected 
by sieving bottom substrate of several types of 
freshwater habitats, mostly lakes, rivers, ponds, 
and wetlands (Figs. 8—11) and preserved in 80% 
ethanol. For comparison, additional sphaeriid 
specimens were included from 39 localities in 
parts of Asia (especially locations in the south of 
the Himalayas, six species from Nepal), Europe 
and Africa (Appendix 1). 

Species determination was performed based 
on shell morphology (soft tissue removed by 
boiling in 3% sodium hydroxide solution) and 
soft body anatomy. Vouchers (DNA, shells) 
are stored at the Systematics and Biodiversity 
collection of the University of Giessen (UGSB; 


Appendix 1). For morphologically determined 
sphaeriids from some Tibetan Plateau loca- 
tions, no molecular data could be obtained 
(Appendix 1). 

Given the conflicting taxonomic concepts 
among Sphaeriidae (Graf, 2013), we here use 
the traditional nomenclature on species level 
(according to Korniushin & Glaubrecht, 2002), 
while classification of the clades inferred is 
based on the phylogenetic analyses of Lee & O 
Foighil (2003). Acknowledging that the purpose 
of this paper is not primarily taxonomic, future 
studies based on more comprehensive global 
datasets should address the internal classifica- 
tion of Sphaeriidae. 


DNA Isolation, PCR, and Sequencing 


This study includes new genetic data for a 
total of 101 specimens, including four speci- 
mens of Eupera spp. as outgroup (Appendix 
1). Genomic DNA was extracted from individual 


FIGS. 8-11. Representatives of the four dominant habitat types inhabited by sphaeriid species on the 
Tibetan Plateau. FIG. 8: Wetland (location CN05); FIG. 9: Lake (CN21); FIG. 10: River (CN26); FIG. 
11: Pond (CNO1). 
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clams using the CTAB protocol described by 
Wilke et al. (2006). For phylogenetic analy- 
ses, we obtained sequences from fragments 
of three genes: (1) mitochondrial cytochrome 
oxidase c subunit | (COI), (2) mitochondrial 
large ribosomal subunit (mtL SU rRNA or 16S) 
and (3) nuclear large ribosomal subunit (nLSU 
rRNA or 28S). Standard and modified prim- 
ers were used for the amplification of COI: 
LCO1490 (Folmer et al., 1994) or a modified 
version (lacking 5 bases from the 5' end) and 
a modified version of HCO2198 (Folmer et al., 
1994; minus the first 6 bases at the 5' end). The 
latter has already been used before especially 
for freshwater bivalves (Hoeh et al., 2002). 
Forward and reverse primers for 16S were 
taken from Palumbi et al. (1991; 16Sar-L and 
16Sbr-H) and for 28S from Park & Ó Foighil 
(2000; D23F and D6R). 

PCR cycling conditions were as follows: an 
initial denaturation step at 95°C for 60 s, fol- 
lowed by 30—36 amplification cycles (denatur- 
ation at 95°C for 30 s, annealing at 52°C for 30 
s, and elongation at 72°C for 30 s), depending 
on the concentration of genomic DNA. The 
reaction was terminated after a final extension 
step at 72°C for 180 s. Sequences (forward 
and reverse) were visualized on an ABI 3730 
XL sequencer (Life Technologies) using a Big 
Dye Terminator Kit (Life Technologies). 

In addition, 16S sequences from 49 sphaeriid 
specimens were obtained from GenBank (Ap- 
pendix 2). 


Sequence Alignment 


The alignments of sequences were done 
individually for all three fragments (68 COI, 
146 16S, and 32 28S). The first base pairs (bp) 
behind the 3’ end of each primer were difficult 
to read in all fragments. Therefore, we trimmed 
each sequence, leaving a 599 bp-long COI 
fragment, a 741—743 bp-long 28S fragment, 
and a 472—503 bp-long 16S fragment. The 
COI and 28S sequences could be unambigu- 
ously aligned using ClustalW implemented in 
BioEdit 7.0.8.0 (Hall, 1999). The alignment of 
16S sequences was performed following the 
instructions for structural alignment by Kjer et 
al. (2009) based on the LSU rRNA structure 
model of Mytilus edulis (Hoffmann et al., 1992). 
In total, three regions were excluded from the 
16S dataset because no reliable alignment 
could be achieved (positions 66—80, 281—296, 
and 341—354), leaving a final alignment of 471 
bp (the aligned dataset is available from the 
corresponding author upon request). 


Phylogenetic Analyses 


The subsequent analyses (except network 
analyses) were performed using all 119 unique 
haplotypes of a combined dataset of COI, 16S, 
and 28S sequences. 

The COI, 16S, and 28S partitions were 
separately tested for substitution saturation 
using DAMBE 5.2.73 (Xia et al., 2003; Xia & 
Lemey, 2009). Under the assumption of a sym- 
metrical as well as an asymmetrical tree, the 
test showed only little saturation for all three 
partitions. 

Phylogenetic relationships were recon- 
structed using Bayesian inference implemented 
in the software package MrBayes v3.1.2 
(Huelsenbeck & Ronquist, 2001; parameters: 
nchains = 4, ngen = 8,000,000, samplefreq 
= 20, temp = 0.1). Prior to the analyses, 
jModelTest 0.1.1 (Posada, 2008) was used 
to find the best fit model of sequence evolu- 
tion for the individual datasets based on the 
corrected Akaike Information Criterion (AlCc). 
For the COI partition the program selected the 
HKY+I+G model and for the 28S partition the 
TIM3+I+G model. Two different models were 
applied for 16S based on the structural align- 
ment: the 16B model (Schóniger & Haeseler, 
1994) for the stem regions and the HKY+G 
model for the loop regions. 

Monitored in Tracer v1.5.0 (Drummond & 
Rambaut, 2007), the combined set of trees 
showed both high ESS (effective sample size) 
values (» 300 for all parameters) and a smooth 
frequency plot. For the consensus tree the first 
25% of the sampled trees (N = 100,000) were 
discarded as burn-in. 


Phylogeographical Analyses 


We performed a parsimony network analysis 
for the most widespread and most abundant 
species P. stewarti using TCS 1.21 (Clement 
et al., 2000) and a combined dataset of COI 
and 16S sequences. The resulting dataset 
contained 14 specimens from nine locations 
across the Tibetan Plateau. 


RESULTS 
Diversity, Distribution and Ecology 
A total of six sphaeriid species were sampled 
across the 27 localities on the Tibetan Plateau 


including five species of Pisidium — P. stewarti 
Preston, 1909; P. subtruncatum Malm, 1855; 


236 


P. zugmayeri Weber, 1910; P. nitidum Jenyns, 
1832; Pisidium sp. | — and one species of Mus- 
culium, M. kashmirensis (Prashad, 1937) (Figs. 
1—6). Pisidium sp. | represents an undetermined 
but morphologically distinct species from the 
clade Cyclocalyx. The Tibetan Plateau taxa 
occurred in different habitats (Table 1). Only 
specimens of P. stewarti and M. kashmirensis 
inhabited all major habitat types: wetlands, lakes, 
rivers and ponds. Lakes were the most species- 
rich habitat type (five species). The most fre- 
quent species on the Tibetan Plateau was P. 
stewarti, found at 21 out of 27 locations (77.8%). 
With a frequency of 29.6% (found at eight loca- 
tions), M. kashmirensis was the second most 
common species. The remaining species were 
relatively rare: P. nitidum and P. zugmayeri 
were each observed at four different locations 
(14.8%); P. subtruncatum and Pisidium sp. | 
were collected at three locations (11.596). At 12 
sampling sites (Table 1) two or more sphaeriid 
species occurred in sympatry (up to four at the 
same location). Atotal of seven different species 
assemblages were found. The most frequent 
composition was M. kashmirensis together with 
P. stewarti (in 5 out of 12 locations; Fig. 7). 
General distribution data of Tibetan Plateau 
sphaeriids were retrieved from the literature 
and the present sampling (Table 2). Pisidium 
nitidum and P. subtruncatum are both known to 
have a large distribution range inhabiting several 
parts of the Holarctic; P. zugmayeri is probably 
endemic to the Tibetan Plateau. The Recent 
distribution of P. stewarti is restricted to areas 
with cold climates in central Asia and Siberia. 


CLEWING ET AL. 


Musculium kashmirensis is distributed in Kash- 
mir, Nepal and on the Tibetan Plateau. Note that 
another species, M. goshaitanensis Nesemann & 
Sharman, 2005, was described from Punyamata, 
Kathmandu, Bagmati, Nepal. 


Phylogenetic Analyses 


The Bayesian analysis revealed a tree with 
five major clades (Fig. 12), which correspond 
to the major lineages of Lee & O Foighil (2003): 
Sphaerium, Pisidium s. str., Cyclocalyx, Odh- 
neripisidium and Afropisidium. The topology is 
unresolved for basal as well as for few terminal 
nodes. Four of the major clades are well sup- 
ported with a Bayesian Posterior Probability 
(BPP) 2 0.95. The genus Pisidium appears 
paraphyletic in the phylogeny as the clade 
comprising Pisidium s. str., including speci- 
mens of P. amnicum (O. F. Müller, 1774) and 
P. dubium (Say, 1817), forms the sister clade 
of the Sphaerium clade, whereas Cyclocalyx is 
the sister group to the former two clades. 

The six sphaeriid species of the Tibetan 
Plateau do not form a monophyletic group 
and belong to the clades of Sphaerium, Cy- 
clocalyx and Odhneripisidium. The clades of 
Afropisidium and Pisidium s. str. are not pres- 
ent on the Tibetan Plateau. The phylogenetic 
relationships of the Tibetan Plateau sphaeriids 
to their respective sister groups are described 
as following (Fig. 12): 

(1) The Tibetan Plateau M. kashmirensis 

are not monophyletic but form a well- 
supported but unresolved clade (BPP 


TABLE 1. Sphaeriid species found among sampling sites on the Tibetan Plateau. Ns = number of spe- 
cies per habitat type; N, = number of locations where the respective species was found; F = frequency 
of species; Habitat Types: wetlands, lakes, rivers, ponds ('others' are described in detail in Appendix 
1); CNO1—CN27 = location codes (see also Fig. 7), * = sympatry occurring. ` 
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TABLE 2. Species of the family Sphaeriidae collected on the Tibetan Plateau with recent distribution 
and range type (except for Pisidium sp. l) based on the literature and own data. 


Species Recent distribution Range type 
M. kashmirensis (Prashad, 1937) Tibetan Plateau1, Kashmir2, Nepal! restricted 
P. nitidum Jenyns, 1832 . Holarctic? widespread 
P. stewarti Preston, 1909 central Asia, north of the Himalayas^5$ restricted (to cold climates) 
P. subtruncatum Malm, 1855 Holarctic? widespread 
P. zugmayeri Weber, 1910 Tibetan Plateau7.8 endemic 


ithis study; 2Prashad, 1937; 3Kuiper et al., 1989; 4Preston, 1909; 5Bogan, 2011; 9Bóssneck, unpubl. data; 7Prashad, 1933; 


8Mischke et al., 2010 


1.00) with M. kashmirensis from Nepal, M. 
japonicum (Westerlund, 1883) from Japan 
and M. miyadii Mori, 1993, from Russia 
(Kuril Islands). The relationship of this 
clade to M. lacustre (O. F. Muller, 1774) 
from Germany, Georgia and the U.S.A., is 
likewise unresolved (BPP 0.96). 

(Il) Pisidium sp. | (clade Cyclocalyx) is not 
monophyletic but clusters with P. caser- 
tanum (Poli, 1791) from Germany (BPP 
0.68). The sister of this clade is a speci- 
men of P. atkinsonianum Theobald, 1876 
from Nepal (BPP 0.62). Note also that P. 
casertanum is not monophyletic according 
to the phylogeny of the present study. 

(II) The Tibetan Plateau specimens of P niti- 
dum (clade Cyclocalyx) form no monophyl- 
etic group but cluster together with other P 
nitidum specimens and specimens of three 
other species (P. edlaueri Kuiper, 1960; P. 
milium Held, 1836; and P. pseudospha- 
erium J. Favre, 1927) from Germany, 
Macedonia and Turkey in an unresolved 
clade (BPP 0.52). 
(IV) The Tibetan Plateau P. subtruncatum 
(clade Cyclocalyx) are monophyletic (BPP 
1.00) and sister to other specimens of this 
species from Georgia, Italy and Macedo- 
nia, together forming a well supported 
clade (BPP 0.99). 
Pisidium zugmayeri (clade Cyclocalyx) 
forms a monophyletic group (BPP 1.00) 
that falls into a strongly supported clade 
of P. casertanum from Ireland, Macedonia 
and Turkey (BPP 0.99); thereby, the speci- 
men from Ireland is the direct sister to the 
Tibetan Plateau specimens (BPP 1.00). 
(VI) The specimens of P. stewarti represent a 
distinct monophyletic group (BPP 1.00) 
within the Odhneripisidium clade, being 
sister to P. moitesserianum (Paladilhe, 
1866) (BPP 0.71). 


3 


The genetic distance between P. stewarti and 
its sister group; however, appears to be high as 
revealed by long branches from their common 
ancestor. For P. stewarti a haplotype network 
analysis was conducted (Fig. 13). 


Phylogeographical Analyses 


The parsimony network analysis of 14 P. 
stewarti specimens from the Tibetan Plateau 
revealed nine unique haplotypes for the com- 
bined dataset of 16S and COI. These haplo- 
types fall into two subclades (A: 6 haplotypes, 
B: 3 haplotypes) separated by eight mutational 
steps (Fig. 13). Specimens from four geograph- 
ically very distant locations (max. ca. 1,400 km) 
share the most frequent haplotype (subclade 
A in Fig. 13). There is no obvious correlation 
between genetic distances (mutational steps) 
and geographical distances on the plateau 
(map in Fig. 13). 


DISCUSSION 
Diversity, Distribution and Ecology 


Prior to the present study, little was known 
about sphaeriid distribution in freshwater sys- 
tems throughout the Tibetan Plateau. Besides 
the records of P. zugmayeri in lakes Bangong 
Co (western Tibet; Prashad, 1933) and Donggi 
Cona (eastern Tibet; Mischke et al., 2010), 
there were only few other records of P. stewarti 
from the plateau (Preston, 1909; Dance, 1967). 
This study, however, demonstrates that the 
diversity of the Sphaeriidae from the Tibetan 
Plateau is significantly higher than previously 
known, as revealed by the first records of M. 
kashmirensis, P. subtruncatum, P. nitidum and 
Pisidium sp. |. Additionally, a recent drilling 
core analysis by Mischke et al. (unpubl. data; 
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FIG. 13. Statistical parsimony network of P. stewarti from the Tibetan Plateau. The geographic dis- 
tribution of subgroups (A, B) is shown in the map. Locations including more than one haplotype are 


labeled in different grayscales. 


ca. 16,000 years ago; Luanhaizi Lake, Qinghai 
province, northeastern Tibetan Plateau, China) 
has revealed shell fragments of P. obtusale 
(Lamarck, 1818). Hence, additional sampling 
might reveal even more sphaeriid species on 
the Tibetan Plateau. 

Sphaeriids are able to occupy a great variety 
of habitats, however, the highest species diver- 
sity is observed in freshwater lakes, ponds, and 
small rivers (Korniushin, 2006). Some sphaeriid 
species show clear preferences to particular 
habitats (e.g., P. lillleborgi Clessin, 1886: lakes; 
P. supinum Schmidt, 1851: rivers), others are 
generalists inhabiting a wide range of habitats. 
Thus, such generalists seem to be less sensi- 
tive to changes in environmental conditions 


<|- 


(e.g., oxygen content; Korniushin, 2006). This 
fact might be of advantage for generalist spe- 
cies occurring in extreme habitats (i.e., high 
mountain ranges like the Tibetan Plateau) with 
partial seasonal freezing, increased ultraviolet 
radiation and low oxygen level. Interestingly, 
P. stewarti and M. kashmirensis were present 
in all major habitat types; the remaining spe- 
cies appear to be more restricted to particular 
habitats (Table 1). This habitat occupation 
pattern might be due to an increased spe- 
cialization of species at the high elevations of 
the Tibetan Plateau. For instance, P. nitidum 
and P. subtruncatum generally inhabit lotic 
(flowing) as well as lentic (stagnant) waters 
(Kuiper et al., 1989), but they have been only 


FIG. 12. Bayesian phylogram of Sphaeriidae (comprising the genera Musculium, Sphaerium, and Pisidium) 
based on COI, mtLSU rRNA (16S), and nLSU rRNS (28S) sequences. Bayesian posterior probabilities 
(higher than 0.5) are given next to the respective node. Individuals from the Tibetan Plateau are highlighted 
in gray and consecutively numbered with roman numerals (see also Fig. 1). The scale bar represents 
substitutions per site according to the applied model of sequence evolution. Major clades are marked with 
white bars according to the classification of Lee & O Foighil (2003). Specimens (sequences) exclusively 
from GenBank are marked with asterisks (*), such with additional new COI sequence data are marked 
with circumflex accents (^). Outgroups (Eupera spp.), with the specimen codes: MW01/1, NA01/1, and 
UG01/1-2, are not shown. Specimens sharing a haplotype are provided next to the phylogram except for 
CN01/1, CN02/1, CN03/3, CN04/1, CN06/1, CNO7/1—2, CN08/1—-2, CN11/1, CN12/1, CN13/1, CN17/1, 
CN19/1, and CN20/1, which share the same 16S haplotype with CN15/2 and of which no COI data are 
available. For detailed specimen information see Appendices 1 and 2. 
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recorded from lakes on the plateau. However, 


a sampling artifact cannot be ruled out due to 
the small sample sizes (N = max. 7) for each 
habitat type. Nonetheless, the highest species 
diversity was observed for the habitat type 
“lakes” (see above), which may be caused by 
a higher diversity of microhabitats in lakes, or 
by the larger size and higher stability of lakes 
compared to ponds and wetlands. 

Recent studies on Tibetan Plateau fresh- 
water gastropods have detected three to four 
distinct genetic lineages for Radix (Oheimb et 
al., 2011), two lineages for Gyraulus (Oheimb 
et al., 2013) and only one lineage for Valvata 
(unpubl. data). Thus, the diversity of Tibetan 
Plateau sphaeriids seems to be slightly higher. 
A relatively high species diversity of sphaeriids 
is also reported for other high mountain ranges 
such as the Andes in South America, as well 
as Mt. Kilimanjaro and Mt. Kenya in eastern 
Africa (Bóssneck, 2012). General observa- 
tions of the freshwater mollusk fauna in high 
mountainous regions show that the gastropod- 
dominated fauna at lower elevations is replaced 
by a sphaeriid-dominated fauna with increasing 
elevation (Bóssneck, 2012). In conclusion, 
the Sphaeriidae appear to be one of the few 
freshwater mollusk taxa being well adapted to 
high mountainous regions. 


Biogeography of Tibetan Plateau Sphaeriids 


Given that the Tibetan Plateau sphaeriids do 
not constitute a monophyletic group (Fig. 12), 
independent colonizations of the plateau have 
to be assumed. We found at least six indepen- 
dent colonization events; four within the clade 
Cyclocalyx, one for Odhneripisidium and one 
for the clade Sphaerium, each representing a 
single species lineage. Phylogenetic analyses 
of the Tibetan Plateau freshwater gastropod 
genera Radix and Gyraulus revealed at least 
three and two independent colonization events 
of the plateau by these taxa, respectively 
(Oheimb et al., 2011, 2013). Apparently, there 
are decisive differences between these three 
mollusk taxa resulting in a potentially higher 
colonization rate of sphaeriids. One major 
difference possibly causing the increased colo- 
nization rate is the high dispersal capacity in 
sphaeriids (see Introduction). 

In order to infer the biogeographical relation- 
ships of the Tibetan Plateau sphaeriid species, 
knowledge about the zoogeographical affini- 
ties of the respective sister groups are useful. 
Generally, the sister groups are dominated 
by species from northern Eurasia (Palearctic 


region). Based on the general distribution of 
the respective species (Table 2) and the close 
relationship to Palearctic taxa, P. nitidum, 
P. subtruncatum and P. zugmayeri can be 
considered as Palearctic elements of the Ti- 
betan Plateau fauna. The affinity of Pisidium 
sp. |, however, remains questionable due to 
the weakly supported relationships to their 
potential sister taxa, the Palearctic species P. 
casertanum and the species P. atkinsonianum 
(Fig. 12). 

In the case of M. kashmirensis several speci- 
mens from the Tibetan Plateau as well as from 
Far East Russia and Japan (Palearctic region) 
cluster with one specimen from the Oriental 
region of Nepal. In contrast to Gyraulus spp. 
(Oheimb et al., 2013), M. kashmirensis may 
thus have crossed the Himalayan mountain 
range. However, for M. kashmirensis from 
Nepal we could only generate sequence data 
from the relatively conservative 16S gene but 
not from COI. Furthermore, this species has 
been collected in Nepal only in the highly fre- 
quented Kathmandu valley and thus may have 
been recently introduced there. As the holo- 
type of M. kashmirensis was described from 
Kashmir (Prashad, 1937), a region with both 
Palearctic and Oriental elements (Das, 1966), 
the zoogeographical affinity of Tibetan Plateau 
M. kashmirensis remains ambiguous. 

The sister species of P. stewarti is P. moi- 
tesserianum in our phylogeny, however, the 
relationship is characterized by relatively long 
branches. Although locality information for 
the respective P. moitesserianum specimen 
was not provided with the original GenBank 
data, the distribution of this species is known 
to be restricted to Eurasia (Zettler & Kuiper, 
2002), and it thus most likely has a Palearctic 
origin. 

Additionally, Pleistocene fossils of P. stewarti 
described from several parts in Europe (Briggs 
& Gilbertson, 1980; Lozek, 2001) also suggest 
Palearctic affinities for this species. 

In contrast to other Tibetan Plateau spha- 
eriids (i.e., P. subtruncatum and P. nitidum 
with a Holarctic distribution), P. stewarti often 
occurs on high elevations in central Asia, 
with a remarkably wide range across the 
plateau. Therefore, this species was used for 
a more detailed haplotype network analysis 
(Fig. 13). The results of this analysis suggest 
high dispersal rates throughout the Tibetan 
Plateau. This is illustrated by specimens from 
geographically very distant locations (i.e., ca. 
1,400 km) which share the same haplotype. In 
contrast, the maximum geographical distance 
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between locations where identical haplotypes 
have been found is approx. 450 km and 300 km 
for Gyraulus and Radix, respectively (Oheimb 
et al., 2011, 2013). This observation highlights 
the high efficiency of sphaeriid dispersal. Nev- 
ertheless, the Himalayan mountain range also 
represents a crucial biogeographical barrier 
for sphaeriids, as for other Tibetan Plateau 
freshwater mollusks. This might be exemplified 
by Palearctic affinities of most of the Tibetan 
Plateau sphaeriids. Furthermore, this fauna is 
conspicuously different from the adjacent Ori- 
ental Sphaeriidae fauna from Nepal comprising 
ten different and only one shared species (M. 
kashmirensis) (see also Korniushin & Bóss- 
neck, 2003; Nesemann & Sharma, 2005). 


CONCLUSION 


The present study demonstrates that the 
Tibetan Plateau Sphaeriidae fauna is charac- 
terized by a previously unknown high species 
diversity, which is dominated by Palearctic 
elements. A total of six different species have 
been recorded, four of them for the first time 
on the Tibetan Plateau. Furthermore, our study 
supports the hypothesis that the Himalayan 
mountain range represents a crucial biogeo- 
graphical barrier for freshwater mollusks. 
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Spec. 
Code 


CN08/3 
CN08/4 
CN09/1 
CN09/2 


CN10/1 
CN10/2 
CN10/3 
CN11/1 


CN12/1 


CN13/1 
CN14/1 


CN14/2 


CN15/1 
CN15/2 
CN16/1 


CN16/2 


Taxon Locality 

Musculium kashmirensis China, Tibet, spring complex 
north of Lake Manasarovar 

China, Tibet, spring complex 
north of Lake Manasarovar 


Pisidium stewarti 


Pisidium stewarti China, Tibet, 
pond near Cam Co 
Pisidium stewarti China, Tibet, 


pond near Cam Co 
Musculium kashmirensis China, Tibet, Garing Co 
Musculium kashmirensis China, Tibet, Garing Co 


Pisidium stewarti China, Tibet, Garing Co 
Pisidium stewarti China, Tibet, 

wetlands south of Gerze 
Pisidium stewarti China, Tibet, 


wetlands west of Cakaxiang 
China, Tibet, Dajiamang Co 
China, Qinghai, artificial dam 

region near Donggi Cona 


Pisidium subtruncatum China, Qinghai, artificial dam 
region near Donggi Cona 


China, Qinghai, Donggi Cona 


Pisidium stewarti 
Pisidium sp. | 


Pisidium stewarti 


Pisidium stewarti China, Qinghai, Donggi Cona 
Pisidium stewarti China, Tibet, 

river west of Ge'gyai 
Pisidium stewarti China, Tibet, 


river west of Ge'gyai 


Latitude 


30.81483° 


30.81483° 


32.08300° 


32.08300° 


30.76610° 
30.76610° 
30.76610° 
32.06443° 


32.44139° 


29103267" 
35.35419° 


35.35419° 


35.2886 1° 
35.28861° 
32.45972° 


32.45972° 


Longitude 
81 .56236° 
81 .56236° 
83.66180° 
83.66180° 


85.01230° 
85.01230° 
65.01280° 
83:90559" 


82.43805° 


85.74019° 
98.34894° 


98.34894° 


98270116" 
98770116" 
80:89955" 


80.89935° 


Spec. 
voucher # 


UGSB 8412 
UGSB 10003 
UGSB 8376 

UGSB 8377 


UGSB 8378 
UGSB 8379 
UGSB 12570 
UGSB 8383 


UGSB 8384 


UGSB 8388 
UGSB 9087 


UGSB 9088 


UGSB 9089 
UGSB 9090 
UGSB 10006 


UGSB 10007 


COI Gen- 
Bank acc.# 


KF483368 


KF483369 


KF483370 


KF483371 


KF483372 


KF483373 
KF483374 
KF483375 


KF483376 


16S Gen- 
Bank acc.# 


KF483258 


KF483259 


KF483260 


KF483261 


KF483262 
KF483263 


KF 483264 


KF483265 


KF483266 
KF483269 


KF483270 


KF483267 
KF483268 
KF483271 


KF483272 


| 28S Gen- 
Bank acc. Z 


KF483325 


no sequence information 


KF483326 


KF483327 


KF483328 
KF 483329 


(continues) 
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9vc 


(continued) 


Spec. Spec. COI Gen- 16S Gen- 28S Gen- 


Code Taxon Locality Latitude Longitude voucher # Bankacc.# Bank acc.# Bank acc.# 

CN16/3 Musculium kashmirensis China, Tibet, 32.45972° 80.89935? UGSB 12571 no sequence information 
river west of Ge'gyai 

CN17/1 Pisidium stewarti China, Tibet, 30.47651° 90.20139° UGSB 10987 KF483273 KF483355 
river near Nam Co 

CN18/1 Pisidium zugmayeri China, Tibet, 33.62848° 79.76558° UGSB 10989 KF483377 | KF483275 
wetlands north of Bangong Co 

CN18/2 Pisidium stewarti China, Tibet, 33.62848* 79.76558° UGSB 10461 no sequence information 
wetlands north of Bangong Co 

CN19/1 Pisidium stewarti China, Tibet, 33.66933° 80.46586° UGSB 10990 KF483274 KF483352 
inflow Bangong Co 

CN19/2 Musculium kashmirensis China, Tibet, 33.66933° 80.46586° UGSB 11001 KF483378 KF483276 KF483351 
inflow Bangong Co 

CN19/3 Musculium kashmirensis China, Tibet, 33.66933° 80.46586° UGSB 11002 KF483277 
inflow Bangong Co 

CN20/1 Pisidium stewarti China, Tibet, 33.56812° 79.94763° UGSB 10991 KF483278 KF483353 
dam region near Bangong Co 

CN20/2 Pisidium stewarti China, Tibet, 33.56812° 79.94763° UGSB 10992 KF483379 | KF483279 
dam region near Bangong Co 

CN20/3 Pisidium nitidum China, Tibet, 33.56812° 79.94763° UGSB 10396 no sequence information 
dam region near Bangong Co 

CN21/1 Pisidium subtruncatum China, Tibet, Bangong Co 33.57666° 79.86945° UGSB 10993 KF483381 KF483281 KF483331 

CN21/2 Pisidium subtruncatum China, Tibet, Bangong Co 33.57666° 79.86945° UGSB 10994 KF483380 | KF483280 KF483332 

CN21/3 Pisidium sp. | China, Tibet, Bangong Co 33.57666° 79.86945° UGSB 10995 KF483282 

CN21/4 Pisidium nitidum China, Tibet, Bangong Co 33.57666° 79.86945° UGSB 10999 KF483285 KF483354 

CN21/5 Pisidium nitidum China, Tibet, Bangong Co 33.57666" 79.86945° UGSB 11000 KF483384 ` KF483286 

CN22/1 Pisidium zugmayeri China, Tibet, 39:92711^ 79.77509" UGSB' 10997 1493882 KF483283 KF483333 


wetlands north of Bangong Co 
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Spec. 
Code 


NO01/1 
NP01/1 


NP02/1 
NP03/1 


NP03/2 
NP03/3 
NP03/4 
NP04/1 
NP04/2 


NP05/1 
NP06/1 


TR01/1 
TRO2/1 
TRO2/2 


UG01/1 
UG01/2 


UG02/1 
VN01/1 
ZA01/1 


Taxon 
Pisidium lilljeborgii 
Pisidium kuiperi 
Pisidium atkinsonianum 
Musculium kashmirensis 
Pisidium annandalei 
Pisidium kuiperi 
Musculium indicum 
Pisidium sp. 
Pisidium atkinsonianum 
Pisidium sp. 
Pisidium annandalei 


Pisidium cf. nitidum 
Pisidium casertanum 
Pisidium personatum 


Sphaerium sp. 
Eupera ferruginea 
Eupera ferruginea 
Pisidium nevillianum 


Pisidium costulosum 


Locality 


Norway, Nordland 
Nepal, Karnali 
Nepal, Karnali 
Nepal, Bagmati 
Nepal, Bagmati 
Nepal, Bagmati 
Nepal, Bagmati 
Nepal, Mechi 
Nepal, Mechi 


Nepal, Mechi 
Nepal, Seti 


Turkey, Isparta 
Turkey, Antalya 
Turkey, Antalya 


Uganda, Mukono 
Uganda, Mukono 


Uganda, Buliisa 
Vietnam, Ninh Binh 
South Africa, Mpumalanga 


Latitude 


68.25166° 
29.08333° 


29.30000° 
27 5968835" 
27:59693" 
27.596895" 
2759693" 
2T 51067" 
27516067" 
27.48717° 
29.47972° 
3/.88664^ 
37.10960° 
37.10960° 
0.22240° 
0.22240° 
2.21527" 
20.25961° 


Longitude 


19185990" 
82.65000* 


82.28333° 
85.37914° 
85.37914° 
85.37914^ 
85.37914* 
87.70000* 
87.70000* 
87.74000° 
81.12778° 
30.858572" 
30:58551" 
30:58551* 
32.609455" 
32.63433° 
31.6%599° 
105.72194° 


-27.01971" 29:87442° 


Spec. 
voucher # 


UGSB 10598 
UGSB 10570 


UGSB 10571 
UGSB 10973 


UGSB 9094 
UGSB 9097 
UGSB 9098 
UGSB 10577 
UGSB 10594 


UGSB 10583 
UGSB 7004 


UGSB 9064 
UGSB 9067 
UGSB 9069 


UGSB 4856 
UGSB 4854 


UGSB 5398 
UGSB 10986 
UGSB 10599 


COI Gen- 
Bank acc.# 


KF483406 
KF483407 


KF483408 


KF483409 
KF483410 
KF483411 
KF483412 
KF 483413 
KF483414 
KF483415 


KF483416 


. KF483418 


KF483417 


KF483419 
KF483420 


KF483421 


KF483422 


16S Gen- 
Bank acc.# 


KF483305 
KF 483306 


KF483307 
KF 483308 


KF483309 
KF 483310 
KF483311 
KF483312 
KF483313 


KF483314 
KF483315 


KF483317 
KF483316 
KF483318 


KF483319 
KF483320 


KF483321 
KF483322 
KF483323 


28S Gen- 
Bank acc.# 


KF483343 


KF483330 
KF483344 
KF483345 


KF483346 | 


KF483347 
KF483348 


KF483349 


KF483350 
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APPENDIX 2. List of additional 16S sequences obtained from GenBank including accession 
number, taxon, specimen code, locality information, and respective author(s). 


16S Gen- 
Bank acc.# 


AF462027 
AF152029 
AF152030 
AF152031 
AF152032 
AF152033 
AF152034 
AF152035 
AF152036 
AF152037 
AF152038 
AF152039 
AF152040 
AF152045 
AF152046 
AF152047 


AF152048 
AY093550 
AY093552 
AY093553 
AY093554 
AY093557 
AY093599 
AY093562 
AY093565 
AY093566 
AY093569 
AY093570 
AY093571 
AY093573 
AY792320 
AY957824 
AY957832 
AY957836 
AY957855 
AY957858 
AY957880 
AY957884 


Taxon 


Pisidium dubium 
Pisidium compressum 
Pisidium variabile 
Pisidium adamsi 
Pisidium sterkianum 
Musculium securis 


Musculium argentinum 


Musculium lacustre 


Musculium partumeium 


Sphaerium corneum 


Sphaerium rhomboideum 


Sphaerium fabale 
Sphaerium simile 
Sphaerium striatinum 


Sphaerium occidentale 


Sphaerium 
novaezealandiae 


Sphaerium tasmanicum 


Musculium japonicum 
Musculium lacustre 
Musculium miyadii 


Musculium transversum 


Pisidium casertanum 
Pisidium casertanum 
Pisidium hallae 
Pisidium nipponense 
Pisidium nitidum 
Pisidium supinum 
Pisidium parvum 
Pisidium japonicum 
Sphaerium nucleus 
Sphaerium corneum 
Pisidium fallax 
Pisidium hallae 
Pisidium ventricosum 
Pisidium milium 
Pisidium nitidum 
Pisidium ferrugineum 
Pisidium ferrugineum 


Spec. 
Code 


US01/1* 
US02/1* 
US03/1* 
US04/1* 
ARO1/1* 
US05/1* 
AR01/2* 
DE06/1* 
US06/1* 
DEO7/1* 
US07/1* 
US08/1* 
US09/1* 
US10/1* 
US11/1* 
NZ01/1* 


AUO1/1* 
JP01/* 
US12/1* 
RZ01/1* 
US13/1* 
DE08/1* 
US14/1* 
N/A.1* 
JP02/1* 
N/A.2* 
N/A.3* 
JP01/2* 
JP03/1* 
UAO1/1* 
EROT” 
US15/1* 
US16/1* 
US17/1* 
US18/1* 
US19/1* 
US20/1* 
US21/1* 


Locality 


USA, Michigan 
USA, Michigan 
USA, Michigan 
USA, Michigan 


Argentina, Buenos Aires 


USA, Michigan 


Argentina, Buenos Aires 


Germany 
USA, Michigan 
Germany 
USA, Michigan 
USA, Michigan 
USA, Michigan 
USA, Michigan 
USA, Michigan 
New Zealand 


Author(s) 


Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 
Cooley & Ó Foighil, 2000 


Australia, New South Wales Cooley & Ó Foighil, 2000 
Japan, Ehime Prefecture Lee & Ó Foighil, 2003 


USA, Michigan 
Russia, Kuril Islands 
USA, Michigan 


Germany, Ammerbuch 


USA, Michigan 
N/A 


Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 


Japan, Gunma Prefecture Lee & Ó Foighil, 2003 


N/A 
N/A 


Lee & Ó Foighil, 2003 
Lee & Ó Foighil, 2003 


Japan, Ehime Prefecture Lee & Ó Foighil, 2003 
Japan, Nagano Prefecture Lee & Ó Foighil, 2003 


Ukraine 
Lithuania, Vilnius 
USA, Colorado 
USA 

USA, Colorado 
USA, Colorado 
USA, Colorado 
USA 

USA 


Lee & Ó Foighil, 2003 


Petkeviciute et al., 2006 


Guralnick, 2005 
Guralnick, 2005 
Guralnick, 2005 
Guralnick, 2005 
Guralnick, 2005 
Guralnick, 2005 
Guralnick, 2005 
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(continued) 
16S Gen- - Spec. 
Bank acc.# Taxon Code Locality Author(s) 
DQ062627 Pisidium N/A.4* N/A Steiner, 2006 (unpubl. 
moitesserianum data) 

EU559087 Pisidium amnicum ALO2/1* Albania, Tušemišta, Schultheiss et al., 2008 
Lake Ohrid 

EU559089 Pisidium amnicum MEO1/1* Montenegro, Crnojevica Schultheiss et al., 2008 
River 

EU559090 Pisidium annandalei INO1/1* India, Pithoragarh Schultheiss et al., 2008 

EU559095 Pisidium casertanum MK05/2* Macedonia, Pretor, Schultheiss et al., 2008 


Lake Prespa 
EUB59097 Pisidium casertanum ESO1/1* Spain, Picos de Europa Schultheiss et al., 2008 
EU559104 Pisidium clarkeanum NP07/1* Nepal, Ghatgai w Sauraha Schultheiss et al., 2008 


FU559105 Pisidium clarkeanum NPO08/1* Nepal, Haripur Schultheiss et al., 2008 

EU559110 Pisidium edlaueri MKO9/1* Macedonia, Struga, Schultheiss et al., 2008 
Lake Ohrid 

EUS559111 Pisidium edlaueri MK10/1* Macedonia, Ohrid town, Schultheiss et al., 2008 
Lake Ohrid 

EUB59133 Pisidium nevillianum NP09/1* Nepal, Lamahi Schultheiss et al., 2008 

EU559147 Pisidium DE03/2* Germany, Uhyst Schultheiss et al., 2008 

pseudosphaerium 


EU559164 Pisidium waldeni NOO02/1* Norway, Davggejavrre Schultheiss et al., 2008 


